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ABSTRACT 

A solar active region (AR) is a three-dimensional magnetic structure formed in the convection 
zone, whose property is fundamentally important for determining the coronal structure and 
solar activity when emerged. However, our knowledge on the detailed 3-D structure prior to 
its emergence is rather poor, largely limited by the low cadence and sensitivity of previous 
instruments. Here, using the 45-second high-cadence observations from the Helioseismic and 
Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), we are able for the 
first time to reconstruct a 3-D datacube and infer the detailed subsurface magnetic structure 
of NO A A AR 11158 and to characterize its magnetic connectivity and topology. This task is 
accomplished with the aid of the image-stacking method and advanced 3-D visualization. We 
find that the AR consists of two major bipoles, or four major polarities. Each polarity in 3-D 
shows interesting tree-like structure, i.e. while the root of the polarity appears as a single tree- 
trunk-like tube, the top of the polarity has multiple branches consisting of smaller and thinner 
flux-tubes which connect to the branches of the opposite polarity that is similarly fragmented. 
The roots of the four polarities align well along a straight line, while the top branches are slightly 
non-coplanar. Our observations suggest that an active region, even appearing highly complicated 
on the surface, may originate from a simple straight flux- tube that undergoes both horizontal 
and vertical bifurcation processes during its rise through the convection zone. 

Subject headings: Sun: general — Sun: solar interior — Sun: surface magnetism 



1. Introduction 

The subject of AR magnetic structure in the 
Solar Convection Zone (SCZ, outer ~ 220 Mm 
of the solar radius) is one of the least under- 
stood topics but it is of crucial importance for 
constraining solar dynamo models and explain- 
ing what drives solar activity and space weather. 
It is widely believed that ARs seen on the sur- 
face are magnetic flux-tubes that are being cre- 
ated by the dynamo pr ocess at a depth in the 
SCZ (|Charbonneaull2QQ5[ ) . Subsequently, the flux- 
tubes emerge through the photospheric surface 
giving birth to ARs or sunspots and magnetic loop 
systems in the corona. On the surface, there is a 
high order of regularity on the pattern of AR mag- 



netic polarities, well d escribed by Hale's and Joy's 
laws (IHale et al. 



19191 ). 



There is limited information on the AR three- 
dimensional (3-D) structu re inferred fr om obser- 
vations of AR emergence. IZwaanl ([19870 provided 
a toy-model explanation of AR emergence by at- 
tributing it to the subsurface structure of an Q- 
loop of magnetic flux with a frayed crest that 
breaks through the surface giving rise to the ob- 
served appea r ance of bipo lar ARs. The work by 



Strous et al 



19960 and IStrous fc Zwaanl (pi 
extends this model to include the horizontal dy- 
namics in order to explain the fact tha t many 
flux-tubes emerge in multiple locations. iTanakal 
(|1991l ) studied complex (delta configuration) ARs 



that exhibit rotational proper motions during 
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their evolution (starting from Hale's-and-Joy's- 
laws-incompatible towards being compatible) and 
attributed i t to a knotted 3-D topology. Also, 
iLeka et all (|l996h found evidence that flux-tubes 
emerge kink-deformed by the current they carry, 
thus containing a twisted 3-D structure differ- 
ent from an ^-shaped flux-tube. However, it is 
generally difficult to determine the detailed 3-D 
structure of an AR, due to the limitation of previ- 
ous observations in terms of temporal and spatial 
resolution. 

On the other hand, there has been a consider- 
able amount of theoretical work which have been 
developed over the past four decades trying to at- 
tack the issue computationally (for a review, see 
Fanll2009l ). The models of emergence in the SCZ 



are ( a) the Thin-Flux-Tube model (TFT [Srjruitl 



1981) and (b) the anelastic MHD model ([Goughl 



1969) . While both models work well in the lower 



SCZ, they might not be valid at the top layers 
of the SCZ (that is, 20 - 30 Mm below surface), 
as the flux-tubes are not thin (TFT assumption 
breaks down) and the velocity field is not sub- 
sonic (anelastic approximation breaks down). The 
reason for making such differences lies behind the 
large pr essure gradient at the layers close to the 
surface ([Fanl 120091 : [Steinl l2012h . Thus, theoretical 
works usually split the SCZ into two parts, the 
lower SCZ and the upper SCZ (~ 20 Mm). 

With the improvement of computational power 
and sophisticated algorithms, it has been possi- 
ble for more "realistic" numerical experiments (i.e. 
fully compressible, radiative-convective 3-D MHD 
simulations) to explore the formation of pores and 



sunspot s (e.g. ICheung et al.l l2010b also see re- 
view by Steinl [2012[ ). although the spatial domain 



achieved so far is still very small (the "deepest" 
simulations go down to a depth of 20 Mm). Also, 
radiative-convective MHD models provide no ex- 
planation on large scale characteristics of emerg- 
ing ARs — such as the Joy's law of AR tilts and 
asymmetric foot-point separation — which are re- 
producible by g lobal-scale models, like the TFT 
approximation ([Caligari et al.l fl995l) and anelas- 
tic MHD simulations ([Fan! 120081 ). To this date, 
because of the computational restrictions of our 
current era and the natural complexity of this 
task, there's no global, fully compressible MHD 
model yet capable in probing the evolution of flux- 
emergence throughout the entire SCZ to the sur- 



face (|Fanll2009f ). 

In this letter, we present the implementation of 
an image-stacking technique as a means to recon- 
struct and study the 3-D structure of an emerging 
AR directly from observations and with great de- 
tail. The NOAA AR 11158 presented in Figure [Q 
is known as the one that produced the very first X- 
class flare (X2.2) of the Solar Cycle 24 a nd its ener- 



getics have been studied thoroughly by ISun et al 
(j2012[ ). The photospheric magnetogram images 
show a complex AR as seen from the surface. How- 
ever, with our 3-D reconstruction method, it is 
rather evident that AR 11158 has a much simpler 
origin, i.e. the horizontal and vertical bifurcation 
of a single progenitor flux- tube. This is proba- 
bly the first study of directly reconstructing the 
detailed subsurface 3-D structure and topology of 
solar active regions. 

2. Methodology 

In this study, we used high time-cadence and 
high spatial resolution observations (O'.'5/pixel) 
taken by t he HMI instrument on board the SDO 
spacecraft (|Schou et al.ll2012h . The HMI instru- 
ment is able to take full-disk maps of the line-of- 
sight (LOS) B-field every 45 seconds. The starting 
time, to, of our selected observation period was on 
10-Feb-2011, 00:00:28 UT when the AR 11158 first 
emerged at heliographic coordinates E53°S20° and 
the ending time of the period under study was 
on 16-Feb-2011, 11:18:27 UT when the AR was 
at W30°S20°, well passed the central meridian. 
During this six-day-long period, the AR had gone 
through the emergence phase and fully developed 
into a mature region. For each of the observa- 
tion frames, we performed a geometrical correc- 
tion to get the radial, or normal field B ni from 
the LOS magnetograms. For further processing, 
we selected a cutout of 240 " x 200 " with a guiding 
center following the AR at the Carrington rotation 
rate. Lastly, we rotated and remapped the solar 
sphere from the center of the cut-outs to the solar 
disk center, practically eliminating both the so- 
lar rotation and the projection effects as well, and 
resulting in a good alignment of cut-out images. 
By such pre-processing, we produced a uniform 
dataset ready to be used in our stacking method. 

The image-stacking method works as follows. 
Each of the frames is a map of the B-field at the 



2 



(thin) photospheric layer. We proceed onto mak- 
ing a stack along the time dimension using a 7.5 
minute cadence of the 2-D cutouts. This cadence 
effectively reduces the number of images by a fac- 
tor of 10, from 12 330 images to 1233 images over 
the time period under study. The choice of the 
number of images accommodates the maximum 
computer capacity in both hardware (in partic- 
ular, the memory) and the software. By starting 
with to at the top of the stack (the X and Y di- 
mension) and adding images at later times con- 
secutively at a lower height (the Z dimension), we 
create a 3-D data cube, which can be used to infer 
the 3-D subsurface magnetic structure of the AR 
prior to its emergence. This technique is based on 
the assumption that the subsurface AR emerges as 
a solid body, i.e. the observed flux on the surface 
at each time instance corresponds to one partic- 
ular height of the body. However, we know that 
an emerging AR is subjected to structural changes 
due to turbulence in the SCZ and the near-surface 
processes. In particular, the near-surface pro- 
cesses may dominate the structure of weaker field 
(more in the discussion section). Nevertheless, the 
AR selected in this paper is a strong AR with 
fast flux-emergence, thus making the surface ef- 
fect minimal. As a first-order approximation, the 
velocity of the emergent structure is assumed con- 
stant, thus each frame contributes equally to the 
height of the structure. For better showing the 
emergent structure, we used only ~ 4.4 days worth 
of data (ending on 14-Feb-2011 04:25:17 UT), i.e. 
the first 800 images instead of the full processed 
dataset. The (x, y, z) final dimensions of the dat- 
acube are 480 pix x 400 pix x 800 pix. The dat- 
acube, which is produced in IDL, is imported to 
PARAVIEW, a visualization software package, for 
further processing and inspection. 

3. Results 

3.1. 3-D Topology of the AR 

In Figure [2j we present the resulted reconstruc- 
tion by showing the iso-surfaces in the 3-D data- 
cubes for two representative constant contour val- 
ues (see online material for a video fly-by around 
the two iso-surfaces). One contour level is at 1100 
G, which shows the "skeleton" of the structure, i.e. 
the core structures of each sunspot /magnetic ele- 
ment of the AR, and the other one at 400 G, which 



shows finer structures that envelope the "skeleton" 
structure of 1100 G. Apparently, the AR is com- 
posed of four major magnetic concentrations or 
polarities, as indicated as PI, Nl, P2, N2, respec- 
tively in the bottom panel of Figure [2j The same 
four polarities are marked in panel (c) and (f) in 
Figure [TJ As seen from the 3-D map, the four po- 
larities originate from two bipoles, Pl-Nl and P2- 
N2. The positive polarity PI connects to its corre- 
sponding negative polarity Nl, and P2 connects to 
N2. In other words, the pairs Pl-Nl, P2-N2 form 
two neighboring flux- tube systems. It is easy to 
note that PI and Nl, and P2 and N2 are not clos- 
ing at their apex, since at this location, the field is 
weak and mostly transverse to the LOS. However, 
if we go down to lower B-field iso-surfaces, i.e. 400 
G, we can clearly see an almost closed system of 
adjacent branch- like arches. 

Instead of coherent, solid flux-tubes, we ob- 
serve a very fragmented, branch-like appearance 
in all polarities of the bipoles. On the surface, 
such fragmentation appeared as the continuous 
emergence of individual small magnetic elements. 
However, these small magnetic elements exhibit 
a remarkably ordered, "swarm" -like collective be- 
havior, separating in terms of polarity and coa- 
lescing in 3-D into big "tree-trunks" — i.e. the 
four polarities of the quadrupolar AR. This tree- 
branch-trunk feature signifies a deeper relation of 
all the small magnetic features within a large- 
scale emerging structure. Each polarity consists of 
several branches, which are almost perfectly con- 
nected to the branches in the opposite polarity 
along the flux- tubes. The branches are probably 
caused by a bifurcation process which is further 
discussed below. 

3.2. Bifurcation in Height 

Both bipoles exhibit similar bifurcation along 
the height (or time). An inspection on the 3- 
D data-cubes reveals that, for each bipole, there 
seems to be a two-phase evolution — or, equiv- 
alently, a "grouping" of the individual small 
branches to just two —but larger— groups we dub 
"Mega-Branches" ( "Mega-Branch-a" and "Mega- 
Branch-/^' or for short "MBa" and "MB/3" , with 
"MBa" preceding "MB/3" in time, as illustrated 
in Figure [2j) . A more thorough study of the mag- 
netic topology should include the time-flux profile 
to characterize the temporal evolution in a more 
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quantitative manner as discussed below. 

In Figure[3j we show the magnetic flux-versus- 
time for each individual polarity of the AR, i.e. 
PI, Nl, P2 and N2, and we overplot each polarity's 
unsigned profile in the same graph. According to 
this plot, Pl-Nl (solid lines) is the first one that 
emerges and it is followed by the second bipole, 
P2-N2, about only six hours later (dashed lines). 
From the time-flux profile in Fig.[3j we see that 
for an individual bipole (i.e. Nl-Pl and N2-P2) 
the magnetic flux of its positive and negative po- 
larity is — to a first order — very similar; such 
similarity is in phase throughout the emergence 
period. Also, for both bipoles, two major flux- 
emergence phases/episodes are identified, with the 
initial one being moderate ("a" -Episode) and the 
later one ("/3" -Episode) being a stronger "flux- 
surge" . The major contribution to these two flux- 
emergence episodes is coming from the respective 
adjacent "Mega-branches" of the emerging flux- 
tube, as suggested by the 3-D reconstruction of 
Fig. [21 In each of the bipoles, the Mega-Branch 
that arrives first at the photosphere, "MBa" , ap- 
pears somewhat weaker and fragmented whereas 
the branch that arrives ~ 2 days later, "MB/3", is 
much stronger. Thus, for each bipole we have a 
bifurcation in height, as deduced from Figs. [2j [3l 

In Table 1, we provide the total unsigned flux 
emerged for each episode's polarity along with the 
information on the duration of emergence and the 
flux-emergence rate, selected manually from Fig- 
ure 03 The duration of emergence is the time 
between the onset (annotated green lines) and 
the end-of-emergence (colored asterisks) for each 
episode. The rate of emergence is defined as the 
flux measured for each episode, divided by its re- 
spective duration of emergence. For comparing 
the individual polarities, we present measurements 
for a+/3-Episodes, i.e. unifying the episodes by us- 
ing the onset time of a-Episode's and for the end 
time, the one of /3-Episode's (thus yielding infor- 
mation on the emergence of the individual polari- 
ties). The AR 11158 is a strong AR where emer- 
gence lasted 110 hours with a flux-emergence rate 
of 5.99x 10 16 Mx s -1 , leading to a total emerged 
unsigned flux of 2.4 x 10 22 Mx. Going down to the 
level of the individual polarities, we should be able 
to quantify the similarities seen in the 3-D visu- 
alization. For both bipoles, the /3-Episode emer- 
gence rates are a factor of 3.2 larger as compared 



to the a-Episode rates, on average. Furthermore, 
the onset times for the Episodes of the Bipoles 
N1-P1/N2-P2 are only ~ 6 hours apart, i.e. very 
similar as it can be also seen in the 3-D reconstruc- 
tion of Fig El 

3.3. Bifurcation in Horizontal Direction 

At first sight, AR 11158's magnetic topology 
seems rather complex. However, the tendency for 
collinearity of the four polarities at the bottom 
of the 3-D cube, suggests that both flux-tubes 
might be related, even originating from the same 
parental magnetic flux-tube. This picture is there- 
fore suggesting that a single sub-photospheric flux- 
tube has been bifurcated into two tubes along the 
horizontal direction , as illustrated in Figure HI 
The almost in-phase evolution of the fluxes with 
time of the two bipoles further reinforces this in- 
terpretation. 

Flux-tubes are non-coplanar when they first 
emerge through the quiet sun (see Fig [T] and the 
online videos). However, right before the emer- 
gence of the "flux-surge" in both of the bipoles, 
there is a strong pushing-aside of the already emer- 
gent Mega-Branches ("MBa's") in a manner like 
they seem to "know in advance" about the arrival 
of the stronger "flux-surge" tubes (or "MB/3's"), 
suggesting that the "MB/3's" interact /collide sub- 
photospherically with the trunks of the "MBa" 
tubes. At a later time, i.e. at the bottom of 
the data-cube, the polarities assume the "correct" , 
Joy's law-compliant tilt. 

Further, for each bipole, the asymmetric polar- 
ity separation suggests flux-tubes with an oblique 
A-shape instead of axisymmetric f2-loops. The A- 
shape has its leading leg more tilted away from the 
vertical direction than the trailing leg. This asym- 
metry can be understood in terms of the Coriolis 
force acting on a r ising flux- tube, as discussed by 
ICaligari et al.1 with the TFT approximation 

and also reproduced using anelastic MHD models 
([Abbett et al.ll2QQlh . 

4. Discussion and Conclusion 

In this letter, we presented a novel image- 
stacking technique for reconstructing the 3-D 
structure of buoyant flux-tubes rising through the 
solar surface and forming observed solar active re- 
gions. Sequences of images have been used before 
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to infer the structure of ARs , e.g. iTanakal (jl99lh 
and iLeka et al.l (| 19961 ). However, the previous at- 
tempts are limited to tracking the locations of AR 
centroids with time, but not the entire structure. 
To our best knowledge, this work is probably the 
first true implementation of the image-stacking 
technique to reconstruct the detailed 3-D structure 
of an AR, using advanced visualization software 
and high-cadence high-resolution magnetogram 
data. 

At least in the early stages of emergence, 
the emerging magnetic structures are two non- 
coplanar neighboring bipoles, but a more detailed 
picture reveals a bifurcated structure for both 
bipoles, in the horizontal direction and along the 
height as well. In the low B-field iso-surfaces, mul- 
tiple magnetic arches can be observed to emerge 
in a continuous manner, with the like-polarities 
coalescing with time. The 3-D reconstruction pro- 
vided good evidence that Mega-branches could 
be originating from the same flux-tube below the 
photosphere. Last, we find that there's a dual- 
phase evolution for both bipoles, as suggested by 
both the topology in 3-D and the time-flux pro- 
file of the AR, providing further evidence for a 
bifurcation in height. Observations also indicate 
that the two bipoles have a common origin. The 
two bipoles have a similar topology in 3-D, similar 
temporal evolution in flux-emergence, and most 
significantly, appear almost collinear at the later 
stage of emergence. It is possible that the two 
bipoles are the result of bifurcation of a single 
progenitor flux-tube early in the evolution. 

It is interesting to note that the 3-D topology of 
the AR 11158 —for each bipole as well as overall— 
exhibits all the qualitative characteristics of the 
TFT approximation. The TFT model's successes 
in repro ducing observations qualitatively are well 
known (Calig ari et al. 1995h . The same qualita- 



tive characteristics of the TFT were reproduced in 



incom pressible MHD simulations bv lAbbett et al 



(|2QQll ) by including the Coriolis force due to the 
solar rotation, in order to study its effects on the 
fragmentation of flux-tubes. This simulation re- 
produces the non-axisymmetrical topology that 
arises due to the Coriolis force. From our recon- 
struction it is evident that we have such an asym- 
metry. The fact that we observe it suggests that 
the upper SCZ has not a severe impact in altering 
the magnetic topology of flux-tubes while travers- 



ing the lower SCZ, after they are born. However, 
the small arch-like magnetic "fibers" seen in the 
400 G isosurfaces may be caused by the strong sur- 
face processes due to a large pressure gradient in 
the upper SCZ. 

This study also demonstrates that the image- 
stacking technique is a promising method for 
studying the 3-D structure of ARs prior to their 
emergence. In the future, we will study the mag- 
netic vector 3-D structure by fully using the mag- 
netic vector observations from the SDO/HMI in- 
strument. 
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Fig. 1 — The first six *days of evolution of the AR 11158 as observed with the SDO/HMI LOS mag- 
netograph. The individual polarities are named after which bipole emerged first, e.g. bipole 1, hence we 
name its negative polarity Nl and its positive PI, and with N2 and P2 emerging at a later time. The 
white cross shows the position of the guiding center of the 240" x 200" FOV at a fixed heliographic latitude 
<j) = —20°. Note that the bipoles initially are non-collinear (dashed lines, panel (c)); at a later time they 
become quasi-collinear (panel (f)). 
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Fig. 2. — The 3-D reconstruction of AR 11158 using the image-stacking method on SDO/HMI LOS magne- 
tograms. The flux-tubes shown here are at \B n \ = 400 G (up) and 1100 G (down). On top of each datacube 
is the last HMI LOS frame of the cube, i.e. the bottom frame, on 14-Feb-2011 04:25:57. The presented 
duration of the observation is 100.4 hours. The positive X-axis direction is westward and Y-axis northward. 
The length of the X-axis is roughly comparable with the height of the Solar Convection Zone (SCZ), i.e. 
about 200 Mm. The black loops are grouping the emergence episodes into Mega-Branches-a,/3. 
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Four Polarities Unsigned Flux vs. Time 




11 -Feb 12-Feb 13-Feb 14-Feb 15-Feb 16-Feb 

Start Time (10-Feb-11 00:00:26) 



Fig. 3. — The time evolution of the flux for the bipoles Nl-Pl (solid lines) and N2-P2 (dashed) suggests 
that both have an early emergence phase (i.e. "a") or "front" of flux, followed by a strong flux "surge" ("/?") 
as also can be seen in FigO The Mega-Branches-a,/3 are the major contributors for each episode. Note the 
persistent lagging of the positive, i.e. leading polarities, PI and P2 with respect to the following Nl and 
N2. The green lines denote the onset of emergence episodes. Also, the end of emergence is shown with an 
asterisk (*) in the respective color. 
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Fig. 4. — Model sketch of the emergence process of AR 11158 in the sub-photosphere. The plane is at 
the depth of the bottom of the SCZ that B-fields originally reside. The cylindrical structures shown at 
that plane in the SCZ are the flux-tubes of the toroidal field, generated by the solar dynamo process. For 
flux-tubes created in the South hemisphere during Solar Cycle 24 (like AR 11158), the B-field vector along 
the tubes is directed from West-to-East (here from right-to- left), as dictated by Hale's law of polarity and 
the Babcock-Leighton dynamo theory. The flux-tubes also develop an asymmetric lambda-shape ("A") as 
they rise. 
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Table 1 

Flux-Emergence Rate, Total Flux Emerged and Duration of Emergence 



Bipole Polarity Episode* Rate a Total b Start Time End Time Duration 

a 1.33 1.52 10-Feb 16:27 12-Feb 00:00 31.66 

PI g 3.87 3.41 12-Feb 18:06 13-Feb 18:30 24.50 

1 1.85 4.93 10-Feb 16:27 13-Feb 18:30 74.15 

a 1.57 1.59 10-Feb 15:52 11-Feb 20:03 28.16 

Nl g 4.43 3.27 12-Feb 16:31 13-Feb 13:10 20.50 

1.87 4.66 10-Feb 15:52 13-Feb 13:10 69.28 

a 0.67 0.91 10-Feb 23:24 12-Feb 13:21 37.99 

P2 g 2.64 6.15 12-Feb 13:21 15-Feb 06:00 64.64 

2 1.91 7.06 10-Feb 23:24 15-Feb 06:00 102.63 

a 0.91 1.24 10-Feb 21:40 12-Feb 11:40 38.00 

N2 p 2.94 5.59 12-Feb 11:40 14-Feb 16:30 52.87 

2.09 6.83 10-Feb 21:40 14-Feb 16:30 90.88 

1+2 5.99 23.79 10-Feb 15:52 15-Feb 06:00 110.16 



*The Mega-Branches-a,/3 are the major contributors to the corresponding emergence 
episodes. 

a xl0 16 Mxs- 1 

b xl0 21 Mx 

c hours 
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